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ABSTRACT 

Quasi-periodic propagating disturbances are frequently observed in coronal 
intensity image sequences. These disturbances have historically been interpreted 
as being the signature of slow-mode magnetoacoustic waves propagating into 
the corona. The detailed analysis of Hinode EUV Imaging Spectrometer (EIS) 
timeseries observations of an active region (known to contain propagating dis- 
turbances) shows strongly correlated, quasi-periodic, oscillations in intensity, 
Doppler shift, and line width. No frequency doubling is visible in the latter. 
The enhancements in the moments of the line profile are generally accompanied 
by a faint, quasi-periodically occurring, excess emission at ~100 km/s in the blue 
wing of coronal emission lines. The correspondence of quasi-periodic excess wing 
emission and the moments of the line profile indicates that repetitive high- velocity 
upfiows are responsible for the oscillatory behavior observed. Furthermore, we 
show that the same quasi-periodic upfiows can be directly identified in a simul- 
taneous image sequence obtained by the Hinode X-Ray Telescope (XRT). These 
results are consistent with the recent assertion of De Pontieu & Mcintosh (2010) 
that the wave interpretation of the data is not unique. Indeed, given that several 
instances are seen to propagate along the direction of the EIS slit that clearly 
show in-phase, quasi-periodic variations of intensity, velocity, width (without fre- 
quency doubling), and blue wing enhanced emission this dataset would appear 
to provide a compelling example that upfiows are more likely to be the main 
cause of the quasi-periodicities observed here, as such correspondences are hard 
to reconcile in the wave paradigm. 

Subject headings: Sun: corona — Sun: UV radiation — line: profiles — waves — solar 
wind 
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Introduction 



Intensity oscillations with a period of three to thirty minutes have been frequently 



observed in polar plumes 



199c 



2000 
20061: 



2OO2I: 



(e.g. 



Ofman et al. 1997: DeForest &: Gurman 



1998 



Baneriee et al.lboiO ) a nd active region loops (e . g.jBerghmans &: Clettel 



Marsh et al. 



Robbrecht et al.ll200ll: Marsh et al.ll2003l : iKing et al 



2003 



19991 



Ofman et al. 



De Moortel et al. 



McEwan fc De Moortel 



2009t Istenborg et al.l I2OI1I ) In spectroscopic studies, small-amplitude 



oscillatio ns (usually a few per cent of the background emission) have been found in line in- 
tensities ( jBanerjee et al.ll2009l ) that are often accompanied by sn iall fluctuations (at most 
a couple of km/s) i n the Doppler velocities ( Wang et al. 2009al lbl: iKitagawa et al. 2010; 
Mariska fc MuglachI I2OIOI ). These quasi-periodic disturbances usually show propagating 



speeds of 50-200 km/s and are almost interpreted as slow-mode magneto-acoustic waves 
propagating into the corona along the magnetic field without exception. 



velocities of 50-150 km/s 


are prevalent in active regions ( 


Sakao et al. 


2007; 


Hara et al. 


2008|: 


De Pontieu et al. 


2009; 


Mcintosh & De 


^ontieu 


2009a 


He et al. 


2010; 


Guo et al. 


2010 


Peter 


2010 




Brvans et al. 


2010 


), quiet Sun ( 


VIcIntosh & De Pontieu 


200c 


)b) a 


nd coronal holes 
as weak upward 


(De Pontieu et al. 


2OO9I: 


Mcintosh et al. 


2010a 


b). These upfiows appear 



propagating disturbances in coronal images, and in spectroscopic observations they are iden- 
tified as significant blue-wing asymmetries in emission profiles of spectral lines formed at 
transition-region and coronal temperatures. These faint upfiows, with a life time of 50- 
150 seconds, are believed to be associated with type-II spicules or rapid blue-shifted event s 



observed in the chromosphere (IDe Pontieu et al.ll2009t iRouppe van der Voort et al.ll2009l ). 



are suggested to provide 
1 coronal heating process 


lot plasmas into 1 


he corona and may thus play an important 


(De Pontieu et al. 


2009 


Mcintosh & De Pontieu 


2009b: 


Peter 


De Pontieu & Mcintosh 


2OI0I: 


Hansteen et al. 


2OI0I). 



These ra pid upfiows often recur at the sa me location on time scales of three to fifteen 
minutes (e.g.. iMcIntosh fc De Pontieul l2009al Jbl). an d would naturall y cause quasi-periodic 
low-contrast oscillations in coronal images (also see IXia et al.l 120051 ). Thus, the discovery 
of these rapi d quasi-periodic upfiows challen ges the universal wave interpretation of coronal 
oscillations. iDe Pontieu fc Mclntoshl ( 120101 ) analyzed time series data that were previously 
studied by IWang et al. fl2009bh to illustrate the presence of slow-mode waves, and found 
that the intensity and Doppler shift oscillations are also accompanied by oscillations in the 
line width and excess emission of the blue wing. They concluded that while the "flows vs 
waves" picture was not unambiguously resolved, the presence of these multi-moment in- 
phase oscillatory signatures is consistent with propagating quasi-periodic upfiows causing 
the observed signature. 
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In this paper, we present new results derived from a timeseries data set obtained by the 
EUV Imaging Spectrometer (EIS. ICulhane et al.l 120071 ) onboard Hinode. We show that co- 
herent oscillatory behaviors in intensity, Doppler shift, line width, and blueward asymmetry 
are clearly present almost everywhere at the root of fan-like structures in the boundary of 
an active region. We demonstrate that these oscillation signatures are caused by repetitive 
high-speed upfiows, which can be direct ly identified in the image sequence simultaneously 
obtained by the X-Ray Telescope (XRT. iGolub et al.l 120071 ) onboard Hinode. 



2. Data reduction and analysis 



The EIS sit-and-stare data used here was acquired in AR 10942 from 17:50 to 20:45 on 
2007 February 20. The 1" x 512" slit was used for the observation, with a 30 s exposure and 
32 s cadence. After standard correction and calibration of the EIS data, a running average 
over 3 pixels along the slit was applied to the spectra to improve the signal to noise ratio. We 
selected two strong emission lines in the spectral window for our study: Fe xil 195. 12A and 
Fe XIII 202. 04A. The Fe xil 195. 12A line is known to be blended with the line Fe xil 195. IsA 
line that is a few percent of the 195. 12A brightness and, since it belongs to the same ion, 
it should exhibit the sam e Doppler behavio r at approximately lOOkm/s in the red wing of 



the line obser yed by EIS (Young et al. 



ignor ed (e.g., iHarra et al 



2008 



20091) . It is common practice that t his weak blend is 



Wang et al.ll2009bl : iKitagawa et al.ll2010l : iTian et al.l 12008 



20101) but this may not always be a safe practice - in regions where hot loop material is 
present the line can be significa ntly impacted by the wings of Fe XIV (195.246A) or an 
unidentified line (see. Table 2 of iBrown et al.ll2008l ): therefore great caution must be taken 
when using this line , especially when c onsid ering possible profile asymmetry analysis of the 
type introduced by iDe Pontieu et al.l (120091 ). Again, as is common practice, we applied a 
single Gaussian fit to each EIS spectrum to derive t he line intensity, Doppl er shift and line 
width for each emission line. Using the method of iKitagawa et al.l ( 120101 ). we performed 
a cross-correlation analysis between intensities of two exposures to obtain the shift of the 
spectra, and thus to remove the jitter in the y-direction. In this case the jitter in the 
x-direction is neglected as it is comparable to the slit width. 



Results of the single Gaussian fit and "R-B" (red-blue) line asymmetry analysis (IDe Pontieu et al 



20091 ) of the Fe xil 195. 12A and Fe xill 202. 04A timeseries are shown in Fig. [21 By assuming 
the average Doppler shift of each line is zero over the entire space-time domain, we calcu- 
lated the relative Doppler shift. We also calculated the non-thermal width of each line profile 
under the assumption of ionization equilibrium (ion temperature equals the formation tem- 
perature of the emission line). To quantify the asymmetry of each line profile, we performed 



-4- 




-550 -500 -450 -400 
X (arcsec) 



Fig. 1. — An XRT image of the studied active region at 19:40 showing the fan-hke structure 
at its boundary. The white vertical hne represents the location of the EIS slit. The section 
between the two horizontal bars representing the root of the fan is selected for further 
analysis. The arrow points to a strong isolated upflow event that travels across the EIS slit, 
see, Fig [21 The online edition of the journal has a animation of this figure. 
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a R-B analysis: after performing the initial single Gaussian fit to measure the line centroid, 
we interpolated the line profile to a spectral resolution ten times greater than the original 
one, then we simply subtracted the blue wing emission integrated over a narrow spectral 
range from that at the same position and over the same range in the red wing. The range of 
integration is then sequentially stepped outward from the line center to build an R-B profile. 
The example presented in Fig. |2] shows the average of the R-B asymmetry from 80-140 km/s 
(normalized to the peak intensity at each pixel). Negative and positive values indicate asym- 
metries in the blue and red wings, respectively. Data in the time range between 122 and 127 
minutes are affected by the spacecraft passage through the South Atlantic Anomaly (SAA). 

At the same time, XRT continuously observed the AR in the Ti-Poly filter with a varying 
exposure time from 2 to 16 s. We selected those frames with an exposure time larger than 
10 s for direct identification of upflow events. In total there are 99 images, with a mean 
cadence of 106 s. These images were first coaligned using a cross-correlation technique and 
then interpolated into regular time intervals. The location of the EIS slit in the image was 
determined by cross-correlating the EIS Fe xiv 264. 78A line intensity along the slit and 
the XRT intensity at different x-locations. The accuracy of the coalignment is about 1" in 
the y-direction. Figure [T] shows an XRT image obtained at 19:40. The white vertical line 
represents the EIS slit location. 

A comparison between Figured] and Figure [2] suggests that strong blue shift, large non- 
thermal broadening, and excess blueward emission are clearly present at the root of the 
fan-like structure in the boundary of the AR (>-83"). We selected the section between the 
two dashed lines in Figure [21 where the values of both Doppler shift and R-B are nega- 
tive throughout almost the entire observation duration, for a more detailed analysis. The 
temporal evolution of the XRT intensity and line parameters of the Fe xil 195. 12A and 
Fe XIII 202. 04A lines in this section are presented in Figure [31 Here the parameters are 
de-trended at each slit location, this is done by subtracting a ten minute running average of 
each parameter from the timeseries. The de-trended intensities and R-B are then normalized 
to the local intensity. The values of Doppler shift and R-B are inverted such that a large 
positive value indicates more blue-shifted emission, or a stronger blueward asymmetry. In 
Figure [H we plot the contours of large non-thermal width on the R-B map of the two lines. 
The data between 122 and 127 minutes were influenced by SAA and thus have been removed 
from the plot. 

We also applied a "R-B guided" double Gaussian fit to the EIS s pectra with a high signal 



P 

to no ise ratio and a strong blueward asymmetry (see, e.g.. Sect. 5 of iDe Pontieu fc Mcintosh 



20101). We used the centroid of the R-B asymmetry profile as an initial guess of the spectral 



position of the secondary component. It turns out that a blue-shifted secondary component is 
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Fig. 2. — Temporal evolution of the line peak intensity, Doppler shift, Non-thermal width, 
and R-B asymmetry (normalized to the peak intensity at each pixel) of the Fe xil 195. 12A 
and Fe xill 202. 04A lines. For the calculation of the Doppler shift, the rest wavelength of each 
line was determined by setting the average Doppler shift as zero. For the R-B asymmetry, 
negative and positive values correspond to blueward and redward asymmetries, respectively. 
The arrow points to an isolated strong upflow event. The section between the two dashed 
lines representing the root of the fan is selected for further analysis. 
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Fig. 3. — Temporal evolution of the de-trended XRT intensity and line parameters of 
Fe XII 195. 12A and Fe xill 202. 04A in the section between the two horizontal bars in Fig- 
ure [H The values of Doppler shift and R-B are inverted so that a large value indicates 
more blueshifted emission or stronger blueward asymmetry. All parameters presented are 
de-trended in the way described in the text. 
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Fig. 4. — Contours of large non-thermal width (top 1/3) superposed on the R-B maps of 
the Fe xil 195. 12A (upper) and Fe xill 202. 04A (lower) lines. Darker and lighter blue colors 
represent stronger and weaker blue asymmetries, respectively. 




Fig. 5. — Timeseries for Fe xil 195. 12A and Fe xill 202. 04A, averaged from y=-54" to -50". 
The black, red, green, violet, blue, and cyan curves represent the line intensity, Doppler 
shift, non-thermal width, R-B, intensity ratio between the secondary and primary Gaussian 
components, and the velocity difference of the two components, respectively. Note that the 
values of Doppler shift and R-B are inverted so that a large value indicates more blue-shifted 
emission or stronger blueward asymmetry. For the non-detrended parameters, the Doppler 
shift, non-thermal width, and velocity difference are shown in km/s, while other parameters 
are shown in arbitrary unit. The non-thermal width is offset by -65 on the y-axis, and the 
velocity difference is divided by 2, for the purpose of illustration within one frame. The R-B 
is also shifted on the y-axis and its zero line is drawn at y=-20, so that values above this 
line indicate blueward asymmetries. For the de-trended parameters, the intensity, R-B, and 
intensity ratio are shown in relative amplitude (percentage), while the Doppler shift and 
non-thermal width are shown in km/s. The intensity, R-B, and intensity ratio are offset by 
8, -15, and 22 respectively on the y-axis. 
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clearly present in the region where strong blueward profile asymmetry is found. An example 
of timeseries for both Fe xil 195. 12A and Fe xill 202. 04A is presented in Figure [51 The 
double Gaussian fit was not performed for the Fe xil 195. 12A line since the weak blend at 
the red wing could have an important impact on the result produced by the sensitive fitting 
algorithm. For a better illustration, we have smoothed the parameters over three adjacent 
time steps. Again, the values of Doppler shift and R-B are inverted. Again, the data between 
122 and 127 minutes have been removed from the plot. 



3. Results and Discussion 



Applying a single Gaussian fit to the line profiles, it has been found that co ronal emission 



lines usually show blue shift of t he order of 30 krn / s at boundaries of sorne ARs (IMarsch et al. 



2004 I2OO8I : iHarra et all I2OO8I : IPel Zannal I2OO8I : IPoschek et all I2OO8I : iMurrav et al 



These blue shifts were thought to b e genesis of the slow solar wind (ISakao et al 



2010h . 



2007; 



Harra et al.ll2008l : IPoschek et al.ll2008l ). However, from Figure[2]we can see that the Fe xil 195. 12A 
and Fe xill 202. 04A profiles at the root of the fan-like structure (y=-83"~-43") in the AR 
boundary are actually very asymmetric with p r ominen t excess emission in blue wings, con- 
firming the results of Mcintosh k De Pontieul (j2009al ). Although the Fe xil 195. 12 A line 
is potentially blended, we place great confidence in its strong blue wing asymmetry here 
since the potential blends are sitting in the red wing of the profile. Such a result suggests 
the presence of continuous fast-moving upfiows (around 100 km/s) and that the emission 
consists of multiple components. The centroid of the line profile derived by using a single 
Gaussian fit only refiects the ensemble velocity of the multiple emission components. And 
different components may have different velocities, which naturally broadens the profile. 

The lack of blue asymmetry above y=-43" seems to suggest that rapid upflows occur 
mainly at the root of the fan-like structures. Perhaps the larger field inclination with height 
along the loops (thus larger angle between the fiow direction and the line of sight) makes 
it difficult to resolve the upfiow signatures. While the profiles of Fe xill 202. 04A are noisy 
there, the Fe xil 195. 12A line profiles reveal clear prominent red wing asymmetries, which 
might be caused by the blends of the line (see above). The lower part of the slit (y<-83") is 
dominated by loop structures in the AR core, where weak red wing asymmetries are present 
in profiles of both lines. It might be related to the complexity of the emission and magnetic 
structures. The asymmetries are more prominent for the Fe xil 195. 12A line, which should 
be related to the blends. 



Examining the XRT movie associated with Fig. [H we can see plasma moving outward 
rapidly along loop structures in the fan. These continuous upfiows have previously been 
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identified by ISakao et al.l ( 120071 ) having an average speed of ~100 km/s. We believe tfiat 
tliese rapid upflows are responsible f or the blueward asymme t ry in coronal emission line 
profiles, consistent with the results of iMcIntosh &: De Pontieul ( l2009a[ ). An isolated strong 
upflow event (indicated by the arrow in Figure [Hand Figure [2]) was clearly visible in both the 
imaging and spectroscopic observations. As the rapidly moving plasma crossed the slit at 
around 19:40, we immediately observed a significant enhancement of the emission in the blue 
wing of the emission line. This enhanced blue wing asymmetry resulted in an enhancement 
in the line intensity, Doppler shift, and non-thermal width derived by a single Gaussian fit. 
This we, believe is an isolated example of the process occurring frequently at the roots of 
the fan structure and is highly unlikely to be the result of wave passage. 

Although the cadence of the XRT observation is much lower than that of EIS, and 
the plasma sampled by the XRT Ti-Poly filter has a higher temperature than the formation 
temperature of Fe xil, from Figure [3] we see a good correspondence in the temporal evolutions 
of the XRT and coronal line intensities. Figure [3] shows that the evolutionary patterns 
of the de-trended intensity, Doppler shift, and non-thermal width are highly similar (the 
correlation coefficient between each pair of line parameters is around 0.6). The R-B pattern, 
to some extent, is also similar to those of the other four. In Figure H] we find that patches 
of large non-thermal width (contours) often coincide with those of large blue asymmetry 
(darker colors). The correlation coefficient between non-thermal width and R-B is 0.40 for 
Fe XII 195. 12A and 0.29 for Fe xill 202. 04A. The relatively low value for the latter might 
be partly caused by the much lower signal to noise ratio of the Fe xill 202. 04A line profiles. 
The coherent behaviors revealed in Figure [3] and Figure H] strongly suggest that continuous 
upfiows with quasi-periodic enhancement of the fiow intensity are responsible for the quasi- 
periodic enhancement of the line intensity, Doppler shift, and non-thermal width determined 
from a single Gaussian fit. We note that the correlations at several instances are not obvious 
or even not present. This is likely to be caused by the poor spe ctral resolution and high 
photon noise of the EIS instrument (IDe Pontieu fc Mclntoshll2010l ). 

As an example. Figure [5] shows the timeseries for Fe xil 195. 12A and Fe xill 202. 04A 
at y=-54"~-50". Correlated changes in intensity, Doppler shift, non-thermal line width, and 
R-B are clearly present. It is also clear that the intensity ratio between the secondary and 
primary Gaussian components generally varies with R-B, suggesting that the fast-moving 
plasma is resolved by our guided double Gaussian fit. Note that the less-than-ideal corre- 
lations at some instances are actually the result of the poor spectral resolution and photon 
noise of the EIS instrument (IDe Pontieu fc Mclntoshl 120101 ). The relative velocity of the 
secondary component is rather stable at ~ 100 km/s, except for several instants when the 
R-B values are relatively small. 
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As we mentioned in the introduction, quasi-periodic intensity oscillations have been al- 
most universally interpreted as waves. However, from Figures [SHSl it is cle ar that repetitive 



flows c an also pro duce oscillatory s i gnatu res, as demonstrated previously by lDe Pontieu fc Mcintosh 



( 120101 ) ■ Recently, IVerwichte et al.l (120 lOf ) presented a slow wave model to argue that the wave 
interpretation is still valid for the observed quasi-periodic intensity perturbations. However, 
the Figure 3 of their paper shows a frequency-doubling of the line width oscillation com- 
pared to the intensity and Doppler shift oscillations, which is not observed by EIS in our 
observation. Instead, all of these single Gaussian parameters show a reasonable correlation 
over several hours. Moreover, in the wave scenario the R-B values often can be positive and 
negative over an oscillation period - this is also not the case in our observation. From our 
Figures [2] and [5] we can see that in the fan root region the R-B values almost remain the 
same sign in the entire 175-minute observation period, indicating the presence of continuous 
blueward emission from upflows. 

In conclusion, we flnd that coherent oscillatory behaviors in intensity, Doppler shift, line 
width, and blueward asymmetry are clearly present almost everywhere at the root of the 
fan-like structures in the boundary of an active region. With coordinated imaging observa- 
tion, we conclude that the quasi-periodicities we observed are more likely to be caused by 
quasi-periodic high-speed upflows. There is no doubt that both waves and flows are present 
on the Sun. We emphasize that it is difficult to distinguish between upflows and waves only 
through the intensity evolution. Spectroscopic observations reveal more information, and a 
combination of imaging and spectroscopic observations is critical for the correct interpreta- 
tion. So far we have found two observations where flows seem to be a better interpretation 
for the quasi-periodicity. With more detailed spectroscopic observations we are sure that we 
can flnd more evidences of flows. 



EIS and XRT are instruments onboard Hinode, a Japanese mission developed and 
launched by ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC (UK) as 
international partners. It is operated by these agencies in cooperation with ESA and NSC 
(Norway). Hui Tian is supported by the ASP Postdoctoral Fellowship Program of National 
Center for Atmospheric Research. The work carried out in this manuscript was partly sup- 
ported by NASA (NNX08AL22G, NNX08BA99G) and NSF (ATM-0541567, ATM-0925177). 
The National Center for Atmospheric Research is sponsored by the National Science Foun- 
dation. We thank Alfred de Wijn and the anonymous referee for helpful comments. 
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